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ABSTRACT 

The current hierarchical merging paradigm and ACDM predict that the z ^ 4 — 8 universe should 
be a time in which the most massive galaxies are transitioning from their initial halo assembly to 
the later baryonic evolution seen in star-forming galaxies and quasars. However, no evidence of this 
transition has been found in many high redshift galaxy surveys including CFHTLS, CANDELS and 
SPLASH, the first studies to probe the high-mass end at these redshifts. Indeed, if halo mass to stellar 
mass ratios estimated at lower-redshift continue to z ^ 6 — 8, CANDELS and SPLASH report several 
orders of magnitude more M ~ halos than are possible to have formed by those redshifts, 

implying these massive galaxies formed impossibly early. We consider various systematics in the 
stellar synthesis models used to estimate physical parameters and possible galaxy formation scenarios 
in an effort to reconcile observation with theory. Although known uncertainties can greatly reduce 
the disparity between recent observations and cold dark matter merger simulations, even taking the 
most conservative view of the observations, there remains considerable tension with current theory. 
Subject headings: galaxies: evolution 


1. INTRODUCTION 

Current theory predicts that galaxies begin their ex¬ 
istence as tiny density fluctuations, with overdense re¬ 
gions collapsing into virialized protogalaxies and even¬ 
tually assemble gas and dust into stars and black holes. 
Although the details of these later stages of assembly 
continue to present a difficult theoretical problem, there 
is broad consensus on the mass and redshift distribution 
of halos produce d through an initial collapse and hier¬ 
archical merging (iSheth et al.ll200ll iSoringd et akll^OOSl : 
IVogelsberger et al.1 l2014aH ' For any specific model, in¬ 

cluding the consensus hierarchical merging model with 
w = — 1 dark energy and without exotic dark matter, the 
halo mass function is straightforward to calculate. The 
predicted halo mass function is also potentially a sensi¬ 
tive probe of the effects of dark matter and dark energy, 
with warm dark matte r models inhibiting halo forma¬ 
tion (jLovell et al.ll20lJl and u> > — I dar k energy mod¬ 
els allowing earlier massive halo formation (jO’Hara et al.l 
l200lH[cTadders et al.ll2n07[ l. The consensus in these mod¬ 
els is a rapid evolution in the density of massive haloes 
at z > 4 that should be observationally evident in galaxy 
luminosity and mass functions. 

Until recently, observations of the highest-redshift 
galaxies were limited to a small number o f massive indi¬ 
vidual galaxies and qua s ars at z > 6 ('e.g.. lFontana et id] 
(j2006ll : lMortlock et al.l (|2011ll l. The recent advent of 
high-redshift surveys including the Cosmic Assembly 
Near-infrared Deep Extragalactic Survey (CANDELS; 
iCrogin et al.l (j20IIID and the Spitzer Large Area Sur- 
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vey with Hyper-Suprime-Cam (SPLASH; Capak et al. 
in prep.) now allows us to probe the galaxy luminosity 
and mass functions, and by implication the correspond¬ 
ing halo mass function, across a range of masses from 
z = 4 — 8. In particular, SPLASH has the broad sky 
coverage required to measure the number density of rela¬ 
tively rare, massive galaxies, while CANDELS is a deeper 
survey over a smaller area, more suited to finding lower- 
mass galaxies. 

A range of work has estimated halo masses for 
z > 4 galaxies and pointed out the tension between 
the expected evolution of the halo mass fun ction and 
galaxy l uminosity and ma ss funct ion (e.g., |Lee et al. 
(20Ijl. iHIldebraiMt (|2n09[l . iFinkelstein et al. 

( 20Hi i. [CaDuti et al.l (|2ni5[l L Here we assemble this pre¬ 
vious work and update it with more recent measurements 
of the galaxy luminosity function, our knowledge of the 
z > 4 galaxy population, and stellar synthesis models. 
We begin by showing that the rapid evolution of the high- 
mass end of the halo mass function at z 4 — 8 allows the 
rare, most luminous and/or massive galaxies in large sur¬ 
veys to provide stringent constraints on galaxy evolution 
in the hierarchical merging paradigm. The most straight¬ 
forward prediction comes from the mass and number den¬ 
sity of the most massive halos that have formed as a 
function of redshift. Furthermore, the consensus ACDM 
model combined with our understanding of galaxy evolu¬ 
tion also predicts a sharp drop in the number density of 
high-luminosity galaxies at fixed L and a rapid decline 
in the luminosity of massive galaxies at fixed number 
density towards high redshift. The details of this rapid 
evolution, if observed, will be sensitive to the baryonic 
physics of early star formation and ACDM models. 

In § [21 we show that the most luminous high-redshift 
galaxies apparently lie in halos too massive to have been 
able to have formed via hierarchical merging by those 
redshifts. In §|3]we show the physical properties of high- 
mass galaxies lie on several scaling relations derived at 
lower redshift, so that it is natural to believe they have 
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been inferred correctly. However, because this result re¬ 
lies heavily upon assumptions about the halo mass to 
light ratio, we are motivated to develop alternative tests 
of hierarchical merging based directly on luminosity func¬ 
tions, that allow us to control for systematic errors. Us¬ 
ing these tests, in § |4j we show that high-luminosity 
galaxies do not exhibit the rapid evolution in number 
density at fixed L expected from hierarchical merging. 
We also discuss the expected conversion between lumi¬ 
nosity and halo mass, showing that an implausibly sharp 
evolution in Mnaio/L would be required for consistency 
with theoretical predictions. Several ways of altering 
galaxy formation models in order to produce these “im¬ 
possibly early galaxies” are considered in § [5j We con¬ 
sider the implications of these disagreements in § |6l 

This wor k uses a = (0.704,0. 272,0.728) 

cosmology (iPlanck Collaboration et ^ I2015I1 through¬ 
out. Unless otherwise specified, mass-to-light ratios refer 
to rest-frame UV wavelengths throughout. 

2. THE EARLIEST GALAXIES AND THEIR HALOS 

In the consensus ACDM model, the high-mass end of 
the predicted halo mass function changes rapidly between 
z ^ 8 — 4, with halos containing the most m assive galax¬ 
ies typ ically virializing towards z = 4 (e.g.. iSheth et al.l 
(|2001I) 1. The timespan of 0.9 Gyr over this redshift range 
means that we likely observe these galaxies within at 
most a few dynamical times of their initial assembly. 
Since galaxies are expected to form after their halos as¬ 
semble, the number density of massive systems and its 
redshift evolution can provide a good probe of the ini¬ 
tial formation of the their dark matter halos. The broad 
redshift range over a relatively small amount of time al¬ 
lows for more precise cosmic epoch measurements than 
are easily obtainable at lower redshift. 

However, at these high redshifts there are a limited 
number of tools available for estimating halo masses. The 
most robust are galaxy clustering based methods that in¬ 
fer the ha lo masses from the spatial di stribution o f gala x- 
ies (e.g., iHildebrandt et all ()2009[ 1: iLee et all (j2012ll b 
Clustering methods have the advantage of not requiring 
assumptions about the physical properties of the galax¬ 
ies, but must assume a model for the dark matter from 
simulations. 

Other methods rely on a relationship between luminos- 
ity or stellar masse s estimated from template fitting (cf. 
Illbert et'af] (12013111 to infer halo masse s. “Abundance 
matching” ('e.g.. iFinkelstein et al.l (|2015[ 11 ties a key fea¬ 
ture of a luminosity or mass function, such as the knee of 
the luminosity function, to a feature halo mass function, 
then matches the galaxy density and the dark matter 
halo density to derive halo masses. The main shortcom¬ 
ing of this method is that it will reproduce the halo mass 
function by definition and is sensitive to the point at 
which the two functions are tied together. Finally, one 
can assume a luminosity or stellar mass to dark matter 
mass ratio measured at lower redshifts where better es¬ 
timat ors of dark matter are available (jLeaiithaiid et aP 

nnn. 

In the past, only a handful of z > 4, galaxies were 
known, far too few to estimate the corresponding halo 
mass function using any of these methods. However, 
recent i nfrared observations by the CANDELS c ollab- 
oration (|Duncan et al.l 12014 iGrazian et al.l 1201511 and 



Fig. 1.— Theor etical halo number density as a function of halo 
mass and redshift IjSheth et al.ll2001l: IMurravge et al.1l2013l) for the 
most massive halos at 4 < 2 < 10 (shown as solid lines, with red¬ 
der colors at higher redshift) compared with observational number 
densities of estimated halo masses corresponding to observed star¬ 
forming galaxies at similar redshifts. Halo masses are estimated 
using clustering (triangle), stellar masses converted to halo masses 
using thelow-redshift scaling ratio Mjf/M* ~ 70 (square), and UV 
luminosities converted to halo masses using ratios determined by 
lower-redshift abundance matching (circle), as described in § [21 
for an overall M^/Mq ^ 120Lijy /Lq. These methods all give 
self-consistent number densities that disagree with theoretical ex¬ 
pectations. We discuss these methods and possible sources of error 
using these techniques in more detail in § [S] 

the Spitzer La rge Area Survey with Hyper-Suprime- 
Cam ISPLASHi lSteinhardt et all (|2(114li i provide a much 
larger sample of high-redshift galaxies with stellar mass 
estimates and can now be used to estimate the halo mass 
function at z ^ 6 — 8. 

A subset of these halo mass function measurements are 
compiled in Figure [U using three different techniques de¬ 
pending upon the type of observation used. For points 
derived a clustering analysis (triangles), the halo mass 
function is measured as a direct result of the technique, 
and we report the derived halo mass and number den¬ 
sity. For points derived using photometric SED template 
fitting, the best-fit inferred stellar mass function is con¬ 
vert to halo mass using Mh/M» ^ 70 derived at lower 
redshift (cf. ILeauthaud et al.! (|2012[l b Finally, for UV 
luminosity functions, the monochromatic UV luminos¬ 
ity is converted to stellar mass using the mean stellar 
mass-to-light ratio derived from abundance matching at 
z = 4, where clustering analysis is also available, then to 
a halo mass using Mh/M^ ^70. The overall conversion 
is Mh/Mq ~ VlQLuy/L q. 

A subset of these halo mass function measurements are 
compiled in Figure [2 showing the large and diverging 
disagreement between the theoretical and observational 
evolution of the halo mass function at high redshift. We 
compare them with theore t ical h alo mass functions es- 
timate d using IShe th et al. l (I2001I1 . as provided by HM- 
FCalc (|Murravge et al.l 120131 1. Soecihcallv. we find that 
observational halo mass function estimates correspond 
to a higher number density of massive halos than should 
have been able to form through the rapid collapse and 
evolution of rare, highly-overdense regions (Fig. [!])• 
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One possible explanation is that the Mh/Luv ratio 
in massive galaxies sharply decreases at z > 4, leading 
to overestimated halo masses for high-redshift galaxies. 
If so, we might hope that this rapid evolution should 
be evident from other measured properties of the galaxy 
population. In § [3J we will use scaling relations to probe 
this possibility, instead finding that this rapid evolution 
is not observed. 


3. HIGH-REDSHIFT GALAXIES APPEAR ‘NORMAL’ 

We have shown that there appears to be a sharp 
conflict between the halo mass function inferred from 
observations at z ^ 4 — 8 and the theoretical pro¬ 
duction of halos in the early universe. We term this 
the impossibly early galaxy problem. However, many 
of these estimates critically depend on galaxy spec¬ 
tral energy distribution templates and scaling relations 
derived at lower redshift. As a result, the possi¬ 
bility that these assumptions have bro ken down by 
z 4 — 8 is a major concern (c f. iConrov et al.l 
(1200911: iSchaerer fc de BarrosI (l20l^: iLee et al.l (|2012ll : 
IStark et al.l (|2013ll : ISoeade et al.l (|2014ll l. and such a 
breakdown would provide an easy solution to the im¬ 
possibly early galaxy problem. 

Lower-redshift galaxies are much more robustly under¬ 
stood, and follow several trends and scaling relations, 
both at fixed redshift and in their redshift evolution. 
A good test of whether high-redshift galaxies are truly 
similar to low-redshift galaxies, as has been strongly as¬ 
sumed, is to check whether they lie on the high-redshift 
extrapolation of this behavior. As summarized below, 
we indeed find that early, high-redshift galaxies appear 
to be ‘normal’, including the most massive early galaxies 
that provide the sharpest disagreement with theoretical 
halo production. 

Over the past few decades, there have been a 
series of observational results indicating that many 
processes in the most massive and luminous galax- 


Madau et al.1 119981: 

iPcrez-Conzalez et al.l 20081) and 

ouasar accretion (IRichards et all 120061: ILabita et all 

20091: ISteinhardt & 
less massive galaxi 
forming galaxies is 
fence of a populati 
z ^ 2 — 4 (IClazeb 

ElvisI 120101). occur earlier than in 
es. The early evolution of star 
further supported by the exis- 
on of massive, passive galaxies at 
rook et all 120041: IDaddi et al.l 120051: 

Carollo et al.l 120131: 

Ilbert et al.l 120131: IStraatman et all 

201 4ll. which must 

have been massive star-forming 


galaxies at high redshift si milar to those r eported 
by SPLASH and CANDELS (|Toft et al.l IMll . Fur- 
thermore, Fell/Mg!! ratios in hig h-redshift quasars 
(jBarth et al.1 l200,i : IXurk e t al.l 120 0711 and ages derived 
from template fittinsf dlVIaraston et aLll^lOfl suggest that 
early, massive galaxies indeed form their stars very 
rapidly. 

In contrast, ACDM theory requires that more massive 
haloes generally assemble later than less massive ones. 
So, reconciling the theory and observations requires that 
more massive galaxies must evolve far more rapidly than 
less massive ones. Because the dynamical time increases 
with mass, this rapid evolution would require a change 
in the baryonic processes that dominate star formation 
at lower redshifts. Such a model would predict that the 
redshift z ~ 4 — 8 universe is a period of rapid transition 


between the initial assembly of massive halos and the 
quick growth of their stellar populations. 

3.1. Similar SFR-Stellar Mass Relations 

To date, SPLASH and CANDELS find no apparent 
deviation from properties derived at lower redshift. In¬ 
stead, they find that the trend towards more and more 
massive galaxies evolving earlier continues out to z ^ 
6 — 8. In particular, at 0 < z < 4 it has been shown 
that star-forming galaxies lie on a “main sequence”, 
with a tight correlation between the existing stellar mas s 
and their star-form ation rate (cf. iPeng et all (|2010fl . 
iNoeske et al.l (|2007f) l. An analysis of over two dozen 
studies of star-forming galaxies using different techniques 
for selection, for measuring stellar mass, and for mea¬ 
suring star formation rate shows strong agreement i n 
both slope and redshift evolution (ISoeagle et al.l [201^ . 
Early, high-redshift galaxies lie on the extrapolation of 
this main sequence to z 6 — 7. 

The star-forming main sequence is well-measured for 
lower-mass obje cts at z ^ 5 — 6, wit h agreement be¬ 
tween SPLAS H (ISteinhardt et al.ll2014f) and CANDELS 
(jPuncan et al.ll2014l) . Where stellar masses are available, 
the most massive early galaxies lie directly on the high- 
mass extrapolation of this main sequence. Because the 
stellar mass and star formation rates are measured using 
different wavelengths, systematic errors from incorrect 
template fitting would produce incorrect inferred quan¬ 
tities in different ways, and likely produce a population 
inconsistent with the main sequence. Similarly, arbitrary 
deblending problems would produce outliers lying some¬ 
where arbitrary, and thus likely off this main sequence. 

The redshift evolution of the star-forming main se¬ 
quence is also well understood at 0 < z < 4, and 
the observed 4 < z < 6 star-forming main sequence 
has the properties produced by extrapolating that time- 
dependence. The increase in mass towards high red¬ 
shift of the N most massive star-forming galaxies (one 
component of what has been termed the ’’downsizing” 
problem) is also we l l measured (Iv an Dokkum et a l.ll2010t 

iDuncan et al.ll201^ ISteinhardt fc SDeaglell2014ll . and the 

most massive star-forming galaxies at z ~ 6 lie on the ex¬ 
trapolation of the best-fit time-dependence at 0 < z < 4 

(Fig. m- _ 

At some point stellar templates derived from low red¬ 
shift will cease to be valid, but this has yet to be ob¬ 
served. It might be expected that this occurs only when 
the physics of star formation have changed, perhaps due 
to very low metallicities producing a very top-heavy IMF 
and even Population HI stars. If so, templates may con¬ 
tinue to be valid well above z ^ 6 — 8. At a minimum, it 
seems clear that they have not become catastrophically 
incorrect in the redshift range where impossibly early 
galaxies have been measured. 

3.2. Quasar-Host Galaxy Relations 

A final test of consistency for early, massive star¬ 
forming galaxies comes from quasar accretion. It is 
believed that there is coevolution between a galaxy 
and its central supermassive black hole. Quasar ac¬ 
cretion exhibits similar behavior to star-forming galax¬ 
ies in the s ense that the re _ is “main seq uence” - 

like behavior (|Richards et al.l 120061 : iLabita et al.l 120091 : 
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Fig. 2.— 96 measurements from 32 studies of the star-forming 
main sequence”^ for fixed M* = Mq at 0 < 2 < 6, adjusted 

to lie on a common set o f calib rations following the prescription 
derived by ISoeagle et al.l (120141) using 25 of these studies. The 
horizontal bars indicate the redshift range spanned by each partic¬ 
ular observation, whil e the vertical errors are the true scatter about 
each MS observation dSoeagle et al.lf2014l) . Although these studies 
use many different methods for determining SFR and stellar mass 
(blue = UV, purple = UV-I-IR, red = IR, green = emission lines, 
yellow = SED fitting, black = radio), they all show good agreement 
with a best-fit (log) linear evolution determined at 0 < 2 < 4 (dot¬ 
ted line), and the highest-redshift measurements lie are consistent 
with the extrapolation of that fit out to 2 6. It therefore seems 

likely that techniques for estimating stellar masses and star forma¬ 
tion rates have not become catastrophically incorrect at 2 ~ 6. 

iSteinhardt fc Elvisll201Cll : iShen et ^120111) with the most 
luminous quasars becoming more luminous and having 
larger black hole masses towards higher redshift. Indeed, 
there has historically been an “impossibly early black 
hole” problem similar to the impossib l y earl y gal axies 

f iresen t ed in this work (iMortlock et al.l ^201111: cf. iCyd 
200311 ■ iMadau fc ReesI (|2001ll . iBromm &: LoebI (|200;^ 1. 
although as we discuss in this work the star formation 
problem appears to be more difficult to solve. 

One parameterization of the similarity between the 
redshift evolution of quasars and star-forming galaxies 
is to note that the ratio of M, for most massive star¬ 
forming galaxies to Mbh for the most massive quasars is 
observed to b e approximately 30:1 a t all f ixed redshifts 
(as shown by ISteinhardt fc Soeagld (1201411 from a liter¬ 
ature compilation). The log M^/Mq ~ 11.2 galaxies at 
z 6 in SPLASH have a similar ratio with early, mas¬ 
sive quasars such as the log MBH/Mc^ ^ 10 .08 quasar 
recently reported at z = 6.4 (|Wu et al.ll2015l l. Because 
virial black hole mass estimates have very different sys¬ 
tematic uncertainties than stellar mass estimates, this 
ratio continuing to hold is additional evidence that the 
properties inferred for impossibly early galaxies are likely 
reasonable. 

Certainly at some sufficiently high redshift, these vari¬ 
ous scaling relations and correlations should break down, 
but again this has yet to be observed. The key point 

^ 25 of these studie s were include d in[SBeasle_et_aIJ I1201J ), wit h 
Magnell^^t_^ (|201J),|Heinis_e^^yj20M),|^^^^^^^t_^y_J201Jh 
DuncaiTer'aLr Il2014 ii. I^chreiberet^dl l l2(iT^ [Paim^a et alJ 
(EoH), and lSalmon et ll201Fil l added for this figure. 


is that the observed high-redshift population of star¬ 
forming galaxies follows the expected redshift evolution 
determined observationally from many lower-redshift 
surveys. Even the most massive, highest-redshift objects 
are found and analyzed using the same standard tech¬ 
niques that have been verified to be successful at lower 
redshift, and are remarkable solely because of the dis¬ 
crepancy between observation and ACDM, not because 
of the data. 

4. LUMINOSITY FUNCTIONS AS PROBES OF HALO MASS 

UV luminosity functions are one of the cheapest and 
most robust observations to obtain at high redshift, and 
can be converted to halo masses by assuming a halo mass- 
to-light ratio. Indeed, most of the high-redshift data 
shown in Fig. [T]and all of those at z > 6 ar e derived from 
monochromatic UV luminosity functions (iMcLure et all 
120091 : iLee et aDl2012l : iBouwens et al.ll2015allbf ). the only 
observations available. The main advantage of this tech¬ 
nique over mass functions is that it places most of the po¬ 
tentially large uncertainty into the halo mass-to-light ra¬ 
tio, which can be parameterized in a single model where 
it can be analyzed more clearly. In contrast, mass func¬ 
tions place much of the uncertainty in the details of the 
diverse set of spectral energy distribution models fit to 
each indi vidual galaxy where they can be difficult to un¬ 
derstand (|Conrov et al.ll200'^ . 

In addition to uncertainty in the correspondance be¬ 
tween halo mass and UV luminosity, comparing obser¬ 
vation with theory also requires matching observations 
of galaxies in the midst of star formation with theoret¬ 
ical predictions about the time at which their halo was 
formed. After all, galaxies likely do not emerge fully- 
formed at the moment a halo virializes; the dark mat¬ 
ter merely has to collapse gravitational ly, while baryons 
get a later start at small masses (cf. iHaiman &: LoebI 
(llQQTlU and must clump and cool in order to form stars. 
If this process typically takes, e.g., 300 Myr for every 
galaxy due to dynamical considerations, then a sharp 
rise in the halo mass function as predicted from redshift 
8.0 to 4.0 would lead to a sharp rise in the UV luminos¬ 
ity function for corresponding galaxies from redshift 6.0 
to 3.4. More generally, in a matter-dominated universe, 
dz/dt oc Thus, delays between halo (smaller t) 

and galaxy formation (larger t) lead to the UV luminos¬ 
ity function evolving equally rapidly with time but over 
a smaller range in redshift, producing a more noticeable 
effect observationally (Fig. [3]). 

These two uncertainties - how to match star-forming 
galaxy UV luminosities to halo formation in both mass 
and time - both lie at the heart of the key open ques¬ 
tions concerning early galactic evolution. Determining 
these two parameters requires a better understanding of 
several key aspects of galactic evolution. The time delay 
between assembly and later evolution depends upon the 
relationship between halo formation and the evolution of 
the galaxies that evolve inside those halos. Constraining 
the luminosity to halo mass ratio involves constraining 
the high-redshift initial mass function, dust content, and 
other processes involved in the formation of the first large 
stellar populations made by galaxies. 

Because there is very little degeneracy between shifts 
in mass and time, if both the UV luminosity function 
and halo mass function are well measured independently, 
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Fig. 3. — Number density of IO^^Mq halos (black, solid) from 
theoretical predictions and corresponding number density (i.e., 
point on a luminosity function) for star-forming galaxies (red, 
dashed) for a toy model in which the stars form 300 Myr after halo 
virializtion. There is a characteristic time for forming such halos, 
seen as a sharp rise in number density over time (top) or towards 
higher redshift (bottom). Although the rise in number density of 
halos and (300 Myr later) galaxies is equally fast in time, because of 
the age-redshift relation, the rise in galaxies, taking place at a later 
time and lower redshift, appears sharper with respect to redshift. 
Because the evolution of mass and luminosity functions is typically 
shown in terms of redshift evolution, this means that a sharp rise in 
number density of halos appears to produce an even sharper obse 
rved evolution in the corresponding luminosity function. 

it should be possible to determine both shifts. Indeed, 
such a determination is the core idea behind abundance 
matching, using a combination of theoretical predictions 
and measurements o f number d ensity to ‘match’ halos 
with g alaxies (e.g., ILee et al.l (I2009fl : IBehroozi fc SilkI 
(|20I5ll : lFinkelstein et al.l (1201511 . What we find, however, 
is that this matching process seems to fail at high red¬ 
shift. As shown in Fig. [1] matching halos with galaxies 
would require that galaxies live in halos so small as to be 
unphysical, in extreme cases less than a factor of three 
higher than stellar masses. Further, the sharp high-mass 
evolution in the halo mass function has no observed coun¬ 
terpart in luminosity functions, either at high redshift or 
over a similar period of time at lower redshift. 

4.1. Halo Mass to Light Ratio Evolution 

Ou t t o z = 4.7 , severa l s tudies (iHildebrandt et ahl 
20091: ILee et all l2012t iFinkelstein et al.l 120151 : 

McLure et al.l l20^ lOuchi et al.l 1200911 1 are consis¬ 
tent with a nearly constant halo mass-to-light ratio. 
However, this ratio must evolve at higher redshifts to 
reconcile the UV luminosity function with hierarchical 
merging (Fig. [T]). There are two ways this might occur: 
halos might be more efficient at producing stars than 
previously believed or stars might produce more UV 
light per unit stellar mass than at lower redshift. 

In order to determine the best explanation, consider 
the size of the change in the halo mass to luminosity ratio 
that would be required to reconcile observation with the¬ 
ory. At a fixed monochromatic luminosity of, for exam¬ 
ple, Mi6oo,ab = —21, there is a 1.5 dex decrease in num- 
ber density from z = 4 to 8 (|Bouwens et al.llTolSall . At 


z ^ 4, clustering analysis shows that t his corresponds to 
a hal o mass of log MijaZo/Af© = 12.4 (IHildebrandt et all 
1200911 . In the most (very likely overly) optimistic sce¬ 
nario, in which stars instantly come into existence so that 
there is no time delay between halo and galaxy forma¬ 
tion, the observed number density for Miqqq^ab = —21 
galaxies at z = 8 corresponds to the number density of 
log Mnaio/MQ = 11.6 halos. Thus, a sharp evolution of 
0.8 dex in Muaio/Luv would be required from z = 4 to 
8. We consider in the remainder of this section whether 
this evolution is plausible. 

4.1.1. Stellar Evolution Models 

Evolution in the halo mass to light ratio could be con¬ 
sidered to be dominated by four effects: 

1. Stellar Evolution: In a typical stellar popula¬ 
tion, the luminosity is dominated by massive stars, 
but the stellar mass (which is used to estimate the 
halo mass) is dominated by low-mass stars. As a 
result, the mass to light ratio is higher for older 
stellar populations. Assuming galaxies at z = 8 
have younger stellar populations than z = 4, this 
effect acts in the direction of reconciling observa¬ 
tion with theory. Increasing metallicity associated 
with maturing stellar populations enhances this ef¬ 
fect. 

2. Changing IMF: If the initial mass function (IMF) 
is top-heavy at high redshift, as expected for Pop. 
HI stars, this would again act to decrease the mass 
to light ratio. However, the IMF is expected to be 
similar at all z < 8 ba sed upon seeming ly analagous 
low-redshift systems (|Dias et al.ll2010ll . and there¬ 
fore to have no impact on impossibly early galaxies 

(§[mi). 

3. Evolving Dust Corrections: Very few galaxies 
with photometric z > 4 are selected as highly dusty 
(cf. Illbert etTn (|2010IU . so a strong reduction in 
extinction from current models is not plausible. 
However, if they were dusty, an increase in extinc¬ 
tion would increase the halo mass to monochro¬ 
matic UV luminosity ratio, making the halos con¬ 
taining early galaxies even more massive than cur¬ 
rently measured. 

4. Merging and Time Delays: Clustering and 
merging results in the addition of both halo mass 
and stellar mass to galaxies, as well as additional 
gas that will eventually form stars. A merger of 
two large objects with the same mass-to-light ratio 
will yield a galaxy with that same ratio. However, 
more gradual accretion of both dark matter and 
gas capable of forming stars will produce an imme¬ 
diate increase in halo mass but delayed increase in 
stellar mass, making the typical halo mass-to-light 
ratio larger at z = 8 than at z = 4. This would act 
in the direction of making the problem worse, be¬ 
cause galaxies at the same luminosity would now 
reside in even more massive halos, which have a 
lower number density and late r viria lization time. 
For example, IBehroozi fc SilkI (|2ni5h find that a 
1O^°M0 galaxy at z = 7 should correspond to a 

halo, whereas we have assumed based 
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Fig. 4.— (a) Expected halo mass to monochromatic UV lumi¬ 
nosity ratio, along with the required evolution to reconcile obser¬ 
vation with theory, and (b) resulting corrected halo mass functions 
derived as in Fig. ^with Mhalo/Ljjv evolving due to a stellar 
population starting at low metallicity at 2 : = 12 a nd agi ng along 
the star-forming main sequence, as described in § 14.1.11 Such a 
model would be reasonable given observational constraints, but 
cannot produce agreement beWeen measured UV luminosity func¬ 
tions and simulated halo mass functions. 

upon lower-redshift measurements that its halo was 

below W^Mq. 

The best-case realistic scenario for producing early, lu¬ 
minous galaxies, then, is one in which an isolated stellar 
population is aging as rapidly as possible (#1) and there 
is no evolution in dust (#3) or time delay between in¬ 
fall and star formation (#4). We model the halo mass 
to light ratio (Fig. |4li) from an initial stellar popula¬ 
tion that fo rmed in one rapid burst, perhap s as early 
as z = 12 (iPlanck Collaboration et al.ll2015r ). followed 
by evolution along the star-forming main sequence un¬ 
til observed at z = 4 — 8. If the time dependence of the 
star-forming main sequence at z > 6 can be extrapolated 
fro m the best-fit evolutio n at z < 6 (as determined by, 
e.g.. iSpeagle et all (|2014ll l. it will be nearly constant at 
higher redshifts because of the short gap in elapsed time. 
Thus, SFR oc M°-^, with the stellar age increasing from 
an initial small value to asymptotically approach 50-150 
Myr. depending upon redshift. This is insufficient to 
reconcile observation with theory (Fig. |4 }d). 

The most likely scenarios are even more difficult to rec¬ 
oncile with theory. We model the evolution of the halo 
mass to light ratio (Fig. as above , but further in¬ 
clude a 300 Myr delay Ccf. TWon^ (|2009[1 1 between merg¬ 
ing and resulting star formation motivated by dynamical 
timescales, so that the luminosity function at time t is 
used to predict the halo mass function at t — 300 Myr. 
The halo mass functions in Fig. [SJj are colored to match 
the redshift of the corresponding luminosity function, so 
that the yellow z = 6.00 {t = 0.954 Gyr) observationally- 
derived halo masses are matched with a yellow z = 8.01 
(t = 0.654 Gyr) halo mass function. 

We additionally i nclude dust evolution as in 
iBouwens et al.l (l2015af) from z = 8 to z = 4. This has no 
effect on halo mass estimates derived from clustering and 
from template fitting that already includes extinction. 


Fig. 5.— Halo mass and UV luminosity functions for a population 
evolving due to (1) a stellar population starting at low metallic¬ 
ity at z = 12 a nd agi ng along the star-forming main sequence, 
as described in § 14.1.11 (2) observed evolution in mean extinction 
as a function of redshift and luminosity; and (3) a 300 Myr time 
delay between dark matter accretion and star formation within ha¬ 
los, so that observed UV luminosity functions are compared with 
merger simulation results at higher redshifts. (a) Expected halo 
mass to monochromatic UV luminosity ratio at two different ab¬ 
solute magnitudes along with the required evolution in the halo 
mass to luminosity ratio to reconcile observation with theory. Ob¬ 
served dust corrections are magnitude-dependent, and therefore so 
is the mass-to-light ratio, (b) Resulting corrected observational 
halo mass functions derived as in Fig. [T] with theoretical halo 
mass functions additional shifted by a 300 Myr time delay, so that 
the yellow z = 6.00 (i = 0.954 Gyr) observationally derived halo 
masses are matched with a yellow z = 8.01 (i = 0.654 Gyr) halo 
mass function. Such a model fits well with current observational 
constraints, but cannot reconcile observation with theory. 


but will alter halo masses estimated from UV luminosi¬ 
ties. For the purposes of this model, we take the z = 4 
abundance matching-derived halo mass to UV luminos¬ 
ity ratios as correct, then apply an additional correction 
at higher redshift to account for evolution in the mean 
extinction at that luminosity between that redshift and 
z = 4. The resulting comparison between observation 
and theory (Fig. [S)d) is likely our best model given cur¬ 
rent data, and produces a larger disagreement than the 
models previously discussed. 

Perhaps the most realistic idea along these lines would 
be to create all of the stars directly at z = 8, in which case 
the age-zero stellar population decreases the halo mass- 
to-light ratio sufficiently to produce a match between the 
z = 8 halo mass and UV luminosity function. How¬ 
ever, this would require these galaxies by z = 4 to have 
stellar populations nearly 1 Gyr old (and thus appear 
passive), and there should be a significant population 
of high-redshift starbursting galaxies. Instead, massive 
galaxies at these redshifts are almost exclusively star¬ 
forming (lllbert et alJl20I(i iLee et ani20I2t [llbert et ^ 
I20I3I IStefanon et al.l 1^131 ISteinhardt et al.ll2014D . tvn- 
ically with < 100 Myr old stellar populations. Further¬ 
more, there is an apparent drop in the densit y of extreme 
starb urst galaxies towards higher redshift (|Casev et al.l 
nnn), contradicting the required evolution. 

Even discarding these apparently contradictory obser¬ 
vations, this sharp shift is not quite enough to reconcile 
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observed luminosities with halo mass functions because 
these stellar populations must then age. It is possible in 
this unphysical model to drive the inferred stellar masses 
low enough to avoid a problem directly at the redshift 
when all of the star formation is assumed to take place 
(z = 8 in this example), but it will reappear a short time 
later and is as sharp of a conflict as before by z = 6. 

Another approach is to allow an inconsistency between 
halo masses and infer red baryon masses in ord er to match 
halo mass with light. iFinkelstein et ahl (j2015f l used mass- 
to-light ratios to determine M* and abundance match¬ 
ing to determine Mnaio from CANDELS, which requires 
that Mnaio/Mbaryon appears to evolve from ~ 40 : 1 at 
z ^ 8 to ~ 100 : 1 by z ~ 4, similar to the overall ratio o f 
70:1 by low redshift (iPlanck Collaboration et ahl 1201511 . 
In the absence of exotic high-energy physics, dark matter 
would collapse at least as early as baryons. Therefore, 
it is difficult to explain how this ratio can increase over 
time. 

However, a stellar population aging due to evolution 
along the star-forming main sequence provides an alter¬ 
native explanation. Consider a low-metallicity galaxy 
beginning with a rapid starburst at z = 10 followed by 
main sequence star formation. As the galaxy grows, the 
average age of the stellar population increases as well, 
and M^,/Luv increases from z = 8 to z = 4. Using the 
z = 4 value of M^,/Lijv at all redshifts will correctly es¬ 
timate M* at z = 4, but overestimate M* by a factor 
of 2-3 at z = 8. Assuming that Mbaryon/M^, remains 
constant, correcting the stellar masses would produce 
Mnaio/Mbaryon ^ 100 : 1 at z = 8. Thus, a reasonable 
evolution in the age of stellar populations could recon¬ 
cile abundance matching with our expectation of con¬ 
stant MHaio/Mbaryon and Mbaryon/M^ if the problems 
with relative density and assembly time of the halos are 
ignored. 

4.1.2. Other Explanations 

Are there other ways to rapidly change the halo mass 
to UV luminosity ratio in order to explain the observed 
mismatch? One possibility is to allow for rapid evolution 
in the initial mass function (IMF). The IMF is expected 
to be substantially the same across cosmic time, even at 
Z ~ 0.01, similar to what is seen in the lowest- metallicity 
stellar populations in near by dwarf galaxies (|Dias et al.l 
120101: iFagotto et al.l [To^ . Thus, a much top-heavier 
IMF would either require a new understanding of early 
star formation or even the possibility that at z = 6, there 
remain residual massive stars from early stellar popula¬ 
tions that formed at even lower-metallicities. Because 
massive main sequence stars have very short lifetimes, 
this is also unlikely. 

The other possibility is that the halo mass to stellar 
mass ratio might have evolved. The standard ratio comes 
from a combination of exp ecting that 10% of ba ryons 
have condensed into stars (|Leauthaud et al.ll2012H and 
that there is a 6:1 dark matter to baryonic matter ratio 
([Planck Collaboration et ^I2015fl . To change this ratio 
by 0.8 dex would either require a complete absence of 
dark matter at z = 8 or that nearly 100% of baryons end 
up in stars at high redshift, both of which would likely 
require new physics that alters halo mass functions as 
well. 

In conclusion, a possible solution for the disagreement 


between hierarchical merging and observation is a change 
in our theory of early star formation, so that we cannot 
easily convert between observed luminosities and halo 
mass. Such a solution would be intriguing in its own 
right, with implications discussed further in §|6l 

4 . 1 . 3 . Dangers of Abundance Matching 

Abundance matching compels agreement between ob¬ 
servation and theory, even though it often comes at the 
expense of having a strong physical motivation for such 
a model. As a result, it can be difficult to understand 
exactly what it means for abundance matching to have 
succeeded or failed. Since any pair of continuous func¬ 
tions can be matched in a way that is empirically correct 
at that redshift, in some sense abundance matching will 
always be successful, and it is difficult to select only the 
physically meaningful matches from a technique that in¬ 
cludes no underlying physics. 

Perhaps, then, abundance matching should really be 
thought of as an extension of the physical model that 
produces a halo mass function. Combining that model 
with observations produces additional constraints that 
essentially become part of the theory, and must be pre¬ 
dicted for that theory to be complete. For example, in 
this work we have shown that ACDM requires rapid, ma¬ 
jor changes in the properties of star-forming galaxies be¬ 
tween z = 4 and z = 8. A logical conclusion is that either 
a mechanism for those changes must be incorporated into 
ACDM or the model must be rejected. As shown earlier 
in this section, producing such a mechanism appears dif¬ 
ficult, but there is also a large space of possible models 
that could be developed. 

One of the reasons for that large parameter space is 
that the galaxy luminosity function is a degenerate com¬ 
bination of many different properties. This is why we can 
conclude that ACDM requires a sharp break in typical 
galactic properties between z = 4 and z = 8, but cannot 
specify precisely which properties must change rapidly. 
Many parameters, including the initial mass function, 
stellar population age, star formation efficiency, extinc¬ 
tion and differential clustering between baryons and dark 
matter combine to produce this discontinuous behavior. 
As a result, it is possible to pick nearly any subset of 
these parameters and abundance match in a way that 
avoids rapid evolution, but always at the cost of sharp 
changes in some of the othe r s. 

For example, iTrac et all (120151) report that the star 
formation efficiency is shar ply variable as a functio n of 
both mass and redshift, and IFinkelstein et al.l (120151) in- 
s tead e x press it as an increasing baryon fraction (see also 
() 14. 1.11) . iBehroozi &: Silld (120151) choose to report an over¬ 
all halo mass-to-light ratio, allowing that to vary sharply 
with both mass and r edshift. In an alternative approach, 
iMashian et al.l (I2015D match observation with theory in 
order to make prediction at z > 10 by removing the con¬ 
tinuity requirement, so that galaxies have no consistent 
history, allowing mass to be both added and subtracted 
in any quantities necessary to match halo mass functions 
at each redshift. 

In this work, galactic evolution is expressed with re¬ 
spect to redshift, but because observations at different 
redshifts are also viewing different ranges in both stel¬ 
lar and (presumably) halo mass, we note that in many 
studies this has instead been expressed with respect to a 
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Halo Occ. 

Baryon Frac 

SF Time 

^halo/t^Q 

^ form 

10% 

100% 

Instant 

5 X lO^i 

5.5 

100% 

30% 

Instant 

2 X 1012 

5.5 

100% 

100% 

150 Myr 

5 X lOii 

6.4 

100% 

10% 

(-1.1 Gyr) 

TABLE 1 

5 X 1012 

3.0 


Various combinations of parameters producing the 

OBSERVED NUMBER DENSITY OF 10^'^MqGALAXIES AT 2 = 5.5, AS 
DESCRIBED IN 

combination of both mass and redshift. Similarly, where 
necessary we compare different quantities by assuming 
that galactic scaling ratios determined at lower redshift 
continue to hold, even though sharp changes in those 
relations might also be consistent with all existing obser¬ 
vations. Ultimately, the proper way to express this effect 
depends upon its cause, and as we have shown, the most 
likely astrophysical effects seem incapable of producing 
such rapid evolution given our current understanding. 

5. PRODUCING MASSIVE GALAXIES IN EARLY HALOS 

Having considered the extent to which observations 
of high-redshift galaxies might allow multiple interpre¬ 
tations, it is important to do the same for our theoreti¬ 
cal understanding of hierarchical merging. Although the 
baryonic physics involved in star formation is quite com¬ 
plex and there are multiple definitions of halo mass used 
in describing the results of simulations, there is broad 
consensus on how dark matter behaves and on the num¬ 
ber density and size of the massive halos that they form. 
Explaining an observed number density of galaxies in 
terms of the density of formed halos depends upon three 
parameters: (1) The fraction of halos containing a galaxy, 
or halo oc cupation rate which is m easured to be ^40% 
at z ~ 5 (jHildebrandt et all 12009( 1: (2) The fraction of 
baryons converted into stars, w hich can be parameter - 
ized and is 10% at low redshift (iLeauthaud et al.ll2012h : 
and (3) The amount of time required after virialization 
for those stars to have formed which will translate into a 
Halo-mass to light ratio and is parameterized by stellar 
population models. We display in Table 1 various combi¬ 
nations of these parameters could produce the observed 
number density of 2 x 10 “® Mpc“^ for M* = 
galaxies at z = 5.5. 

We find that each of these combinations requires im¬ 
plausible physics, such as 100 % of baryons being turned 
into stars instantly upon halo virialization or 10 % of 
baryons forming stars over 1 Gyr before the dark mat¬ 
ter halo virializes, or contradicts other observational re¬ 
sults. For example, several combinations require halo 
masses below IO^ ^ M^t^. However, at 3.1 < z < 4.7, 
iHildebrandt et all (12009( 1 used clustering measurements 
in CFHTLS to find that galaxies at 25.5 mag are found 
in halos of log Mhalo/M q = 12.3. A similar ratio would 
yield \ogMh„.in/M(h = 1 2.8 for massive galaxies at z = 6 . 
iFinkelstein et ah] (|2015ll use CANDELS observations to 
argue that the dark matter to baryon ratio decreases to¬ 
wards higher redshift, so that the halo mass to stellar 
mass ratio is 50:1 at z = 6 and 40:1 at z = 7. This would 
correspond to log Mhaio/MQ = 12.7. 

In summary, solutions that give a plausible halo oc¬ 
cupation fraction require an implausibly short timescale 
for star formation, a much higher fraction of baryons to 
be converted into stars than the 10 % in current models. 


or both. Any solution with a standard ratio between the 
halo mass and stellar mass of between 50:1 and 100:1 
{Mhaio = 5 X 10^^ — 10^^Mq) requires most of the star 
formation to occur well in advance of initial collapse and 
virialization. 

5.1. Massive Galaxies in Merger Simulations 

Extensive effort has been put into studying the for¬ 
mation of massive galaxies and their dark matter halos 
through numerical simulations. The vast difference in 
both scale and dominant physical processes between hi¬ 
erarchical merging and star formation means that sim¬ 
ulations cannot investigate both processes directly, but 
rather use semi-analytical prescriptions for connecting 
the properties of massive galaxies to their halos. 

These prescriptions attempt to model baryonic physics 
on a very macroscopic level and are drawn from lower- 
redshift relationships observed between galaxies and 
their host halos. However, extrapolating these rela¬ 
tionships often lead s to unphysical resu lts: the Mil¬ 
lennium simulation (iSprineel et al.l l2005h can produce 
M* = 10 ^^M 0 galaxies at z = 6 , but they live in dark 
matter halos with M = Mq. This is a very small 

halo mass to stellar mass ratio on several fronts: Theoret¬ 
ically, it would require the baryons to cluster in advance 
of much of the dark matter, as well as very nearly all 
baryons to have ended up in stars by z = 6 . Observation- 
ally, i t wo uld be in conflict with t he IHildebrandt et ahl 
()2(I09tl and IFinkelstein et ahl (j2015l l halo mass to stellar 
mass ratios discussed in the previous section. 

The Illustris simulation (iVoeelsberg er et ahl l2014bl : 
iSoarre et ahl 120151 : I Wellons et al.l 1201511 picks a set of 
baryonic relationships that avoids these unphysical ex¬ 
trapolations, resulting in a stellar mass function and 
luminosity function that look similar to the halo mass 
function. As a result, the simulated number densities 
of massive galaxies are consistent with observation out 
to M, ^ 1 O^°M 0 at z ^ 4 — 6 , but very few galaxies 
are produced with M, > 1 O^°'®M 0 , disagreeing with the 
observed number densities at that mass and redshift. 

The conclusion from these simulations is that there is 
broad consensus on the halo mass function, but consid¬ 
erable freedom in matching those halo masses to galactic 
properties. Even given that freedom, matching observa¬ 
tions would require that either stellar masses are vastly 
overestimated or there is a sharp disconnect between the 
baryons in massive galaxies and their dark matter ha¬ 
los. Specifically, simulations are consistent with typical 
galaxies at these redshifts, but cannot produce the earli¬ 
est, most massive galaxies seen at z > 4 with the intro¬ 
duction of reasonable physics. Rather, then, understand¬ 
ing these galaxies appears to require additional physics 
not yet included in these simulations, whether baryonic 
physics relating to star formation or high-energy physics 
altering the timing of massive halo formation. 

6. DISGUSSION 

We have shown that recent observations of high- 
redshift galaxies are inconsistent with current theoret¬ 
ical models of galactic assembly. As a general principle, 
when theory and observation disagree, it is historically 
best to believe the observational result. However, in this 
case the observations also rely on untested theoretical 
assumptions about stellar evolution. Thus, something is 
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wrong, but what? We can divide the possible flaws and 
explanations into three possible categories: 

1. Failed Template Fitting or Redshift Deter¬ 
mination: If the measurements are wrong, it is 
most likely not in the redshifts (as spectroscopic 
redshifts exist for many, albeit less massive galax¬ 
ies at 2 > 5) but rather in the assumption that 
templates derived from lower-redshift galaxies can 
be used at z = 6. This implies that the halo mass 
to monochromatic luminosity ratio changes sharply 
above z = 4. As discussed in § 14.11 the most likely 
explanation for this evolution would be a sharply 
top-heavier IMF at higher redshifts. This can po¬ 
tentially be tested in the near future using super¬ 
nova rates, and certainly following the launch of 
JWST. 

Later spectroscopy confirmed that low-redshift 
templates yielded correct results for inferred quan¬ 
tities such as stellar mass out to z < 3. If 
this breaks down by z 6, it would mean that 
purely photometric surveys are now insufficient, as 
new models must be developed and observation- 
ally tested by JWST. It is inevitable that at some 
point astronomers must encounter this problem, 
but would be unpleasant to discover that it hap¬ 
pens at a lower redshift than currently believed. 

2. New Clustering Physics: Another possibility 
is that halos indeed collapse earlier than allowed 
by current models, something that would simul¬ 
taneously solve the high-redshift massive quasar 
problem as well. Current collapse times are de¬ 
rived from gravity acting on perturbations that can 
be confirmed using cosmic microwave background 
measurements, and therefore a much more rapid 
collapse of the halos would require the introduction 
of new high-energy physics. Possible solutions here 
might then provide exciting new constraints on the 
nature of dark energy or dark matter. It should be 
noted that many dark matter models under current 
consideration are warm dark matter, which would 
suppress the z 6 halo mass func tion rather than 
enhance it (|Gao &: TheunsI 1200711 . Dark energy 
with w > —1 could enhance early structure forma¬ 
tion, although cosmic microwave background ob¬ 
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servations (iPlanck Collaboration et ^ l20I5fl cre¬ 
ate considerable tension with the w > —0.95 re¬ 
quired to solve this problem (|0’Hara et al.l[2005 
[^adders et al.ll2007ll 

3. Early Star Formation: A baryonic solution is in¬ 
stead to allow main sequence star formation much 
earlier than the initial collapse of halos. Such ideas 
would need to evade difficult constraints from both 
low-redshift observations, as well as solve the prob¬ 
lem of cooling to form small stars at low metal- 
licities. We further note that if initial stars form 
in small clumps rather than in fully-formed pro- 
togalaxies, the “clu mpyness factor” in reionization 
(|Ouchi et al.l 1200^ will be much higher than ex¬ 
pected, resulting in far more rapid reionization 
from star formation alone than currently expected. 
If so, this will be evident in infrared background 
fluctuation measurements. 

All three answers carry major consequences for both 
our current understanding of the initial stages of galactic 
formation and our future plans for studying high-redshift 
galaxies. So, better observations are needed. There is 
considerable hope that followup observations will help to 
determine whether the first of these three explanations 
is the right one. 

Rather than speculate as to which explanation is best, 
we instead stress that the high-mass objects coming out 
of high-redshift surveys are now critically important. Fu¬ 
ture surveys should concentrate on finding and charac¬ 
terizing these objects in sufficient numbers to constrain 
how these galaxies and their halos co-evolve. Since the 
earliest, most massive galaxies are rare, wide area surveys 
on the > 1 degree scale will be needed. These objects are 
no longer merely an extra point or two in the last panel 
of a figure, but rather pose a key problem at the heart 
of high-redshift extragalactic astronomy, and need to be 
given corresponding attention. 
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